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利用 帕 囚 托 非 支 配 关 系 实 现 
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摘 要 : 在 多 目标 进化 算法 中 ， 时 间 复 杂 度 过 高 是 普遍 的 


多 运算 时 间 。 和 针对 三 目标 问题 ， 
在 一 定 宛 余 操 作 ， 需 对 全 部 的 解 先 排 序 后 ， 才 能 再 分 配 等 


利用 帕 累 托 非 支配 关系 结合 差分 进化 ， 实 现 高 效 三 目标 进化 算法 。 算 法 每 次 选 代 对 种 群 中 
选择 后 代 个 体操 作 ， 当 后 代 种 群生 成 时 便 跳 出 计算 ， 从 而 减少 个 体 的 计算 数量 ， 
时 给 出 该 方法 的 相关 理论 分 析 和 证 明 过 程 。 然后, 针对 一 系列 


在 分 配 等 级 同时 进行 


一 :一 
[| 


[S] 


效 三 目标 差分 进化 的 方法 


忠义 上-， 盛 梦 园 :， 王 林 景 : 
读 息 工程 学 院 ， 郑 州 450052; 3. 香港 科技 大 学 计算 


p! 


问题 ， 特 别 是 三 个 目标 函数 以 上 时 ， 解 的 等 级 分 配 占 用 了 过 


利用 帕 累 托 支配 关系 ， 对 解 的 等 级 分 配 进行 研究 ， 发 现 经 典 的 等 级 排序 及 分 配方 法 存 


级 并 选择 下 一 代 ， 造 成 部 分 不 必要 的 运算 。 为 减少 该 宛 余 ， 
级 的 个 体 进行 计算 ， 
降低 运算 量 ， 同 
目标 优化 问题 , 将 提出 方法 与 著名 排序 方法 NSGAII， 


已 言 迄 


最 高 等 


及 近年 来 优秀 的 ENS 方法 进行 对 比 实验 。 仿真 实验 结果 表明 ,提出 方法 在 时 间 复 杂 度 和 收 伊 速 度 上 优 于 经 典 方法 ， 稍 
差 于 ENS 方法 。 在 标准 测试 函数 DTLZ1-DTLZ6 的 性 能 上 ， 提 出 方法 近似 于 ENS 方法 ， 优 于 NSGAII 算法 ， 从 而 验 


证 了 提出 方法 的 有 效 性 和 正确 性 。 
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Using Pareto dominance relationship to realize efficient tripe-objective 


differential evolution algorithm 


Xu Yulong!»;, Pan Xu!, Wang Zhongyif, Sheng Mengyuanl, Wang Lingjing! 
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Hong Kong University of Science & Technology Hong Kong, China) 


Abstract: In the Multi-objective evolutionary algorithm, the solution of the rank distribution takes up most of the computation 


time. According to the Pareto dominance relationship, we research the solution of the rank and distribution in this paper, and we 


found that there are some redundancy operation of rank and allocation in the classic methods. The classic algorithms need to 


rank all solutions, and select the next generation after sort, 


so it caused some unnecessary sort operations. To reduce the 


redundancy, we introduce a fast method of non-dominated sorting. The proposed algorithm only deals the individual which 


belong the highest level in current population, and meanwhile it can choose individual into the next generation during the process 


of distribution level. When the next generation individuals were selected enough, the program is ended. This method has reduced 


the number of the sorting of individuals, and the time complexity is also reduced greatly. Then, we use three objectives 


optimization problem as an example, combine the fast non-dominated sorting method with the differential evolution algorithm, 


and propose a fast three objective differential evolution algorithm based on non-dominated sorting. Finally, we compare our 


method with the famous ranking method and ENS/BNS method. The simulation results show that our proposed algorithm in 


time complexity and convergence speed is better than the classical method, and is not worse than the ENS/BNS algorithm. On 


the performance of the standard test function DTLZ1-DTLZ.6, 


the proposed method is similar to or better than the comparison 


algorithm. Based on above experiment, we have verified the correctness and effectiveness of our proposed method. 


Key Words: tripe-objective; Pareto; non-dominated solution sorted; convergence speed 
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多 目 


标 优化 问题 (multi-objective optimization problems， 


标 函 数 一 般 多 于 两 个 , 其 解 是 一 个 集合 , 称 为 非 支 配 


你 


(nondominated seb) 或 Pareto 最 优 解 集 (Pareto-optimal set) 
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。 有 多 种 方法 可 以 求解 一 般 优化 问题 ， 如 线性 规划 、 梯 度 下 
箱 、 进 化 算法 等 ， 但 相关 成 果 表 明 ， 多 目标 进化 算法 (multi- 
objective evolutionary algorithms, MOEAS) 较 适宜 求解 多 目标 优 
化 问题 0 。 


多 


化 受到 国际 学 者 的 广泛 关注 ， 大 量 优秀 的 算法 被 


标 进 
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| ChinaXiv 合 作 期 刊 
录用 稿 许 玉龙 ， 等 : 利用 帕 累 托 非 支配 关系 实现 高 效 三 目标 差分 进化 的 万 法 
提出 ， 其 中 以 NSGAII 算法 外 和 SPEAII 算法 中 最 为 车 名。 上述 ” 类 和 数据 分 析 的 方法 对 多 目标 差分 进化 算法 进行 改进 、 文 献 
两 种 算法 都 是 通过 空间 密度 和 帕 累 托 支 配 来 判断 解 的 好 坏 ， 通 [24,25] 提 出 基于 网 格 索引 指导 进化 方向 的 多 目标 差分 进化 算法 ， 
过 构造 帕 累 托 候选 解 集 来 获得 最 优 解 。 国 内 学 者 中 ， 焦 李 成 、 都 取得 了 较 好 的 效果 。 但 是 ， 上 述 多 目标 差分 进化 方法 在 对 非 
陈 火 旺 、 崔 逊 学 、 郑 建国 、 袭 自 兴 等 人 中 3， 对 多 目标 进化 的 主 。 支配 排序 时 ， 需 要 对 全 部 的 个 体 分 配 等 级 ， 会 造成 部 分 匈 余 操 
要 技术 、 理 论 分 析 及 发 展 历史 进行 深入 研究 ， 并 指明 了 研究 方 作 。 为 改进 该 问题 ， 本 文 前 期 对 两 目标 差分 进化 算法 做 了 相关 
向 和 面临 的 问题 。 其 中 ， 研 究 提高 算法 计算 效率 是 一 个 非常 有 ”的 研究 29I， 介绍 了 一 种 快速 的 非 支 配 解 排序 方法 ， 并 得 到 了 较 
意义 的 问题 。 好 的 效果 ， 但 该 方法 只 是 对 两 目标 优化 问题 的 研究 。 

为 提高 算法 的 计算 效率 ， 学 者 们 进行 了 不 同 的 尝试 ， 此 ， 本 文 在 已 有 工作 的 基础 上 ， 介 绍 利用 帕 累 托 非 支配 
NSGAI 算法 外 相 比 上 一 代 版 本 及 同时 期 其 他 著名 算法 , 在 计算 关系 实现 高 效 标 差分 进化 的 方法 ， 针 对 标 优化 问题 ， 
速度 上 有 较 大 优势 ， 但 其 全 局 时 间 复 杂 度 仍 为 O(n”)。 该 方法 能 够 在 进行 非 支 配 解 排序 同时 ， 选 择 个 体 进入 下 一 代 。 

Tang 等 人 中 依据 帕 累 托 支 配 关系 对 该 问题 进行 改进 ， 使 用 ”在 实验 对 比 测试 中 ， 该 方法 与 经 典 和 最 新 算法 进行 对 比 研 究 ， 
擂台 原则 对 得 到 的 解 分 配 等 级 ， 从 种 群 中 随机 选择 一 个 解 ， 把 里 论 分析 和 实验 结果 都 证 明了 提出 方法 的 正确 性 和 高 效 性 。 
区 en 1 ”相关 理论 
群 中 某 个 解 支 配 “ 擂 主 ”， 用 解 蔡 代 原 来 的 “ 擂 主 ”， 算法 的 时 
间 复 杂 度 和 空间 复杂 度 分 别 为 O(mn?) 和 O(n)。 在 仿真 实验 中 ， ”1.1 差分 进化 算法 
该 方法 展示 出 比 NSGAII 更 高 的 计算 效率 。 同 样 为 提高 计算 效 差分 进化 算法 自 提出 以 来 ,凭借 速度 快 、 易 实现 、 参 数 少 等 
率 ，Clymont 等 人 00 提 出 了 著名 的 非 支配 排序 算法 : deductive 特点， 迅速 成 为 进化 计算 领域 中 最 著名 的 算法 之 一 。 差 分 进化 
sort， 该 方法 能 记录 所 有 解 之 间 的 比较 结果 , 推断 所 有 解 的 支配 遇 于 启发 式 算 法 , 包括 了 变异 、 交 叉 、 选 择 、 更 新 等 基本 操作 ， 
关系 ， 从 而 避免 多 余 的 比较 ， 算 法 的 时 间 复 杂 度 和 空间 复杂 度 ”为 节省 篇 幅 ， 下 面 简要 介绍 其 基本 流程 。 

分 别 小 于 O(mn”) 和 O(n)， 实 验 结果 验证 了 该 方法 的 有 效 性 。 首先 随机 产生 初始 种 群 ， 然 后 依 该 种 群 进行 变异 操作 ， 

不 久 前 ， 张 兴 义 等 人 Hl 提出 一 种 高 效 的 非 支 配 解 排序 方法 ， ”多 种 变异 策略 可 供 使 用 ， 常 见 的 有 
命名 为 ENS。 该 方法 最 大 特点 是 : 当 一 个 候选 解 与 已 分 配 等 级 a) 随机 引导 变异 DE/rand/1: 

的 解 比较 时 ， 得 出 此 候选 解 的 等 级 ， 避 免 了 大 量 匈 余 操 作 。 算 Wi 0) 
法 首先 对 第 一 个 目标 函数 值 按 升 序 排序 ,此 后 ,候选 解 pn 与 F(F | | | 

为 当前 检查 等 级 ) 中 的 解 只 有 两 种 支配 关系 : 候选 解 被 等 级 下 中 b) 最 优 引 导 变 异 DE/ best/1: 

的 解 支配 ; 候选 解 不 被 等 级 下 中 的 解 支配 。 作 者 基于 ENS 提出 站 O) 
了 ENS-SS 和 ENS-BS 两 种 排序 算法 。 在 ENS-SS 中 ， 从 第 一 | . - 

个 等 级 开始 检查 ， 直 到 p, 被 分 配 等 级 结束 ， 如 果 pn 被 五 中 任 c) 当前 到 随机 引导 变异 DE/current-to-best/1: 

意 解 支配 ， 则 继续 检查 后 续 等 级 ; 如果 z 不 被 尺 中 所 有 解 支 人 @) 
配 ， 将 py 分 配 到 五 中 ; 而 在 ENS-BS 中 ， 从 中 间 等 级 开始 ， 检 | 1 

查 范 围 初 始 化 为 已 存在 等 级 。 每 检查 完 一 个 等 级 ， 范 围 缩减 一 其 中 :x 是 目标 向 量 ，* 是 变异 向 量 ， 下 标 i 表示 第 i 个 个 体 ， 
半 。 如 果 疡 被 正中 任意 解 支配 , 检查 范围 缩减 为 F 后 的 所 有 等 。 而 1,2,3,.….,NP,NP 是 种 群 数量 , g 表示 第 g 代 种 群 。 式 3) 中 xpevs 
级 ; 如 果 户 不 被 尺 中 所 有 解 支配 , 检查 范围 缩减 为 之 前 所 有 是 种 群 中 最 优 个 体 ，%sg 是 当代 个 体 ，wsg 是 变异 后 个 体 ，xig, 


等 级 ， 直 到 
算 效 率 上 显著 优 


差分 进化 算法 (differential evolution, DE) 


Pr 被 分 配 到 某 个 等 级 中 。 实 验 表 


于 对 比 算法 。 


出 的 优秀 算法 2 


注 。 文 献 [1 和 较 早 介绍 了 使 月 


一 ， 它 以 高 效 、 


化 问题 ， 
最 著名 的 多 


目标 


取得 较 好 的 效果 。Robic 等 人 提 昌 
杨 平等 人 09 使 用 选举 法 则 来 构 


算法 DEMO 


15] 。 


简便 的 特性 


月 ，ENS 算法 在 计 


“5 是 近年 来 涌现 


受到 了 众多 学 者 关 
昌 DE 结合 非 支 本 


关系 求解 多 目标 


上 了 基于 差分 进化 


造 帕 累 托 非 支 配 解 集 ， 并 进行 理论 证 明 其 能 够 提高 效率 ， 取 得 
了 理想 的 成 果 ， 但 其 成 果 仅 是 对 非 支 配 排序 次 数 的 比较 。 任 长 
安 等 人 0 提出 一 种 改进 的 基于 目标 空间 分 割 的 多 目标 进化 算 
法 ， 算 法 基于 区 间 索 引 值 进行 高 效 的 运算 ， 但 该 方法 仅 在 两 目 


标 问 题 上 进行 了 测试 。 类 似 的 ， 学 者 们 对 差分 进化 算法 及 其 在 


实际 中 的 应 用 也 


做 了 较 多 研究 52 。 不 久 前 ， 文 献 [23] 结 合 聚 


Xizg, Xrag 为 群体 中 随机 选择 三 个 互 不 相同 的 个 体 ，(xig- xizg) 即 


为 差分 向 量 ，A 是 缩放 因子 ， 用 于 对 差分 向 量 缩放 ， 控 制 搜索 
步 长 。 
然后 ， 进 行 交叉 操作 ， 把 生成 的 变异 向 量 vis 与 个 体 向 量 


Xig 进行 交叉 ， 得 到 实验 向 量 wsg。 差 分 进 


化 的 交叉 方式 包括 


项 式 交 又 和 指数 交叉 ， 一 般 情 


如 下 : 
i 
We 
a :8 
WiG = 
好 ， 


其 中 : wi, 表示 G 代 种 群 中 第 i 个 体 的 第 个 基 
E 数 。Cr 6 [0 ,1]， 是 交叉 概率 


是 问题 的 旨 
速度 。 随 机 生成 数 randij 小 于 
从 变异 向 量 继 承 ， 否 则 ， 从 目 


况 ， 


常 | 


if rand, [0,1| < Cr or j= jo 


的 是 二 项 式 交 又 ,方法 


otherwise 人 

姑 ,j=1,…,D,D 

刀子 ， 控 制 算法 收敛 的 

Cr 或 者 j 等 于 jrard 时 ， 实 验 向 量 


标 向 量 继 承 。 Jrand 是 [1,D] 内 的 随 
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china 


录用 稿 


机 数 ， 它 保证 实验 向 量 


变异 向 量 。 


接 下 来 是 选择 操作 , 差分 


方法 是 将 目标 向 旧 


时 与 实验 向 量 


用 高 强度 贪 禁 的 选择 策略 ， 
进行 逐一 比较 ， 优 秀 者 进化 到 下 


一 代 。 以 求解 最 小 化 为 例 ， 其 选择 方法 如 式 (5) 所 示 。 


u, , 让 fw) < f(x) 


Xion = 
i 


以 上 概括 了 差分 进 


献 [12]。 


1.2 多 目标 优化 概述 
在 多 目标 优化 中 ， 各 个 
和 的 ， 问 题 的 最 优 解 通常 是 一 个 解 和 


标 优 化 中 
本 质 的 区 


多 目标 优化 过 


Pareto 支配 关系 ， 


多 个 Pareto 最 优 朋 


最 优 解 是 唯一 的 ， 所 以 多 


程 中 ， 两 个 解 之 间 的 
多 上 
你， 这 些 Pareto 最 


题 的 最 优 解 集 。 解 和 


化 算法 的 简要 流程 ， 其 详 


(5) 


目标 函数 之 间 是 相互 冲突 、 不 可 调 
而 不 是 一 个 解 。 而 在 单 目 


标 优化 与 单 目标 优化 有 


集中 的 所 有 人 解 之 间 都 是 非 支 玫 


占 优 关 系 又 称 为 帕 累 托 
标 优 化 问题 不 存在 全 局 最 优 解 ， 会 存在 
改 合 ， 就 称 为 问 
Cb 关系， 决策 者 


依据 自己 的 爱好 和 需求 ， 从 中 挑选 一 个 或 多 个 作为 最 终 的 最 优 


解 使 用 。 多 目标 进 


化 算 当 


洲 
洁 
ra 
HB 
芒 
3 
Ud 

No 
TH 


长 的 Pareto 最 仇 
数学 描述 参见 文献 [5]，》 


hm 


多 的 、 分 布 均匀 的 、 


标 进化 相关 的 概念 和 
i 空间 在 此 不 再 资 述 。 


2 ”基于 帕 累 托 非 支 配 关 系 实现 高 效 三 目标 差分 进化 


2.1 高 效 的 非 支配 解 排序 方法 


NSGAII 是 最 著名 的 多 目 
都 参考 了 它 的 思想 和 过 程 。 而 


许多 改进 算法 
本 文 发 现 该 算法 存在 计算 匈 余 ， 


算法 中 的 非 支 配 解 排 序 过 程 需 对 所 有 个 体 分 配 等 级 后 ， 排 序 操 


作 才 结束 ， 然 后 才能 进行 选择 
算法 ENS00 中 ， 也 需要 对 所 有 个 体 》 


下 一 代 操 作 。 最 近 被 提出 的 著名 
分 配 等 级 ， 但 ENS 避免 


了 所 有 非 必需 和 部 分 必需 的 支配 比较 次 数 ， 提 高 了 算法 运行 速 


度 。 上 述 算 法 中 ， 都 是 


后 按 个 体 等 级 由 高 到 低 排 序 ， 从 等 级 1 到 等 级 的 顺序 选 出 染 


色 体 进入 下 一 代 


然后 选择 出 N 个 进 
对 部 分 2N 个 体 进行 计算 ， 就 外 
色 体 。 这 将 使 参与 计算 的 个 体 数量 1 
所 提高 。 因 此 ， 依 据 三 目标 


在 多 目标 进化 中 ， 假 设 种 群 
代 的 染色 体 数量 是 即 是 N， 而 经 
算 , 需要 对 所 有 2N 个 染色 体 进 


E 对 种 群 所 有 染色 体 进 


入 子 代 。 设想 在 排 


本 文 提 出 高 效 的 三 


有 2N 个 体 进 行 计算 对 上 


见 算 法 1. 


行 等 级 确定 ， 然 


名 


规模 为 N， 需 要 选择 进入 下 一 
2N 个 染色 体 计 
排序 、 等 级 分 配 、 判断 比较 ， 

序 和 等 级 分 配 时 ， 能 否 仅 
找到 满足 进入 下 一 代 的 N 个 染 
去 的 运行 效率 便 有 
至 排序 的 洪 在 规律 ， 


目标 非 支配 解 排 


序 方法 ， 该 方法 不 需要 对 所 


上 ， 能 够 降低 时 间 复 杂 度 ， 


算法 1 高 效 的 
输入 : 合并 后 的 初始 种 群 和 变异 种 群 U， 


标 非 支配 解 排序 方法 


体 算法 参 


色 体 个 数 为 2 


倍 的 种 群 
输出 : 满足 进入 下 一 代 种 群 的 染色 体 ， 个 数 为 N 
1 for 等 级 rank=1 
2 ”for 个体 二 1:2N 
3 找到 属于 第 rank 级 的 4 个 特殊 个 体 ,存储 到 集合 


4 这 4 个 分 别 是 到 万 轴 、 户 轴 、 方 轴 、 原 点 距离 最 小 


情况 下 ，4 个 个 体 之 间 可 能 会 重复 
对 集合 Srank 中 个 体 分 配 等 级 rank 
/与 集合 Srank 做 差 集 ,得 到 凡 
for 个 体 j=1: 路 
if 个体 j 与 集合 Srank 中 任 一 个 体 都 非 支配 
10 则 j 等 级 为 rank, 将 j 存储 到 Srank 集合 
11 endif 
12 end for 个 体 j=1: 路 
13 得 到 第 rank 等 级 全 部 个 体 ， 存 于 集合 Srank 
14 ”将 集合 Srank 中 个 体 全 部 进入 子 代 
15 放 进入 子 代 个 体 总 数 sum>N 


16 溢出 情况 发 生 , 进行 拥挤 距离 计算 , 选 出 部 分 个 
体 

17 end if 

18 ”进入 子 代 个体 总 数 sum==N 

19 结束 程序 ,得 到 子 代 个 体 

20 end if 


21 end/fori=1:2N 

22 集合 UV 与 Srank 做 差 集 ,结果 记 为 0 集合 

23 令 WV=Q 即 用 0 更 新 

24 ”rankt+, 考 虑 计算 下 一 个 等 级 个 体 

25 end rank 

算法 1 中 ， 令 种 群 规模 W， 算 法 输入 为 2 倍 的 种 群 规模 集 
合 0， 和 输出 为 选 出 YX 个 染色 体 进 入 子 代 。 提 出 的 排序 方法 在 选 
择 下 一 代 个 体 时 ， 无 须 对 所 有 的 2N 个 体 进行 计算 对 比 ， 减 少 
了 个 体 比 较 次 数 ， 避 免 了 存在 的 计算 元 余 。 主 要 特点 如 下 : 首 
先 计 算得 到 当前 种 群 中 ， 属 于 最 高 等 级 的 4 个 特殊 个 体 〈 假 设 
不 重合 )， 将 他 们 存储 到 集合 Srank。 然 后 ， 将 种 群 剩余 个 体 与 
Srank 中 个 体 逐 一 比较 ， 当 个 体 i 被 集合 Srank 中 的 某 个 体 支 配 
时 ， 结 束 与 集合 中 剩余 个 体 的 比较 ， 节 省 了 比较 次 数 。 当 个 体 
i 与 集合 Srank 中 所 有 个 体 都 是 非 支 配 关系 时 ,将 个 体 i 加 入 集 
合 Srank。 在 比较 时 借鉴 文献 [11] 的 经 验 : 从 集合 Srank 中 最 后 
个 体 开始 ， 使 得 每 次 更 新 Srank 后 ， 其 个 体 都 是 按 第 一 个 目标 
函数 值 升序 排列 ， 所 以 Srank 中 靠 后 个 体 更 可 能 支配 个 体 i。 量 
复 该 操作 进行 处 理 , 最终 得 到 当前 种 群 中 最 高 等 级 的 所 有 个 体 ， 
等 级 记 为 ranf 三 1， 存 储 在 集合 Srank。 和 集合 Srank 中 个 体 有 待 
进入 下 一 代 ， 如 果 集 合 中 个 体 全 部 进入 下 一 代 后 ， 下 一 代 个 体 
总 数 超过 NW， 则 进行 拥挤 距离 计算 选择 部 分 个 体 ， 如果 下 一 代 
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孙 用 入 许 玉龙 ， 等 : 和 用 本 条 攻 并 二 到 关 和 DB /全 和 和 同 
个 体 总 数 不 超 过 N， 则 所 有 个 体 进入 下 一 代 ， 并 将 rank 增 1， 2.3 ”高 效 三 目标 差分 进化 算法 框架 和 时 间 复 杂 度 分 析 
更 新 种 群 VU， 重复 上 述 过 程 ， 选 择 出 第 二 等 级 的 个 体 ， 如 果 下 依据 提出 的 高 效 三 目标 非 支配 解 排序 方法 ,结合 差分 进化 ， 
一 代 个 体 总 数 等 于 N， 则 程序 结束 ， 不 需要 对 剩余 的 所 有 个 体 提出 一 种 基于 非 支 配 解 排序 的 高 效 三 目标 差分 进化 算法 


进行 计算 ， 提 高 了 计算 效率 。 根 据 一 般 经 验 ， 到 进化 过 程 的 中 EMODE)。FMODE 的 总 体 流 程 如 算法 2 所 示 。 
期 ， 算 法 只 需 计算 前 几 个 等 级 个 体 ， 便 可 满足 下 一 代 的 选取 ， 算法 2 高 效 三 目标 差分 进化 流程 
从 而 不 需 再 对 剩余 个 体 进行 计算 ， 较 大 地 节省 了 时 间 。 输入 : 算法 参数 、 种 群 规模 N、 进 化 代数 FES 
2.2 高效 三 目标 非 支配 解 排序 算法 理论 基础 输出 : 进化 结束 后 的 解 集 
提出 的 高 效 三 目标 非 支 配 解 排序 算法 中 ， 最 重要 的 是 确定 1 种 群 初始 化 ，X(0) :x7(0),…, xx(0); 
四 个 特殊 的 个 体 ， 并 保证 这 些 个 体 属于 当前 种 群 中 最 高 等 级 。 2 加 载 算法 参数 G, N,Cr ,F. 
需要 说 明 的 是 ,对 于 三 目标 优化 而 言 ,最 多 只 有 4 个 特殊 个 体 ， 3 while( 进 化 开始 不 满足 结束 条 件 ) 
于 这 些 特 殊 个 体 可 能 会 重合 ， 所 以 特殊 个 体 数目 具体 应 为 1 4 ”变异 ， 依 据 式 (1) 得 到 M 
至 4 个， 在 此 本 文 假设 一 般 情况 他 们 不 重合 ， 所 以 需 确定 为 4 交叉 ， 依 据 式 (4) 得 到 O 
个 个 体 。 如 何 证 明 这 些 个 体 在 当前 种 群 中 属于 最 高 等 级 ， 是 较 合并 M 、Q， 得 到 U，U 含 有 2N 个 个 体 


使 用 提出 的 高 效 三 目标 算法 ， 参 见 算法 1 
算法 1 找到 第 rank 等 级 的 个 体 时 , 令 为 天 

10 ”将 KK 个 个 体 设 为 选择 进入 下 一 代 个 体 时 

11 if 己 选择 个 体 总 数 等 于 和 N 

12 break， 进 化 代数 加 1 

13 ”六 已 选择 个 体 总 数 小 于 和 N 

c 14 玉 个 个 体 进入 下 一 代 

0 15 ”已 选 择 个 体 总 数 大 于 和 N 

16 对 第 rank 个 等 级 进行 拥挤 距离 计算 

17 选择 合适 的 个 体 , 进入 下 一 代 , 使 得 下 一 个 个 体 为 


3 
6 
为 关键 的 基础 ， 证 明 过 程 如 下 : 7 ”从 品 中 选择 满足 条 件 的 入 个 个 体 进入 下 一 代 
8 
9 


图 1 高 效 三 目标 非 支配 解 排序 中 四 个 特殊 个 体 


根据 相关 理论 , 给 出 如 下 推理 。 为 描述 方便 , 如 图 1 所 示 ， N 个 
本 文 使 用 x、y、z 分 别 代 蔡 fi/、fp、 有 ff。 18 break， 进 化 代数 加 1 
推理 1 若 存 在 个 体 a(x,y,z) 的 x 上 坐标 值 ， 是 当前 种 群 中 所 19 end while 


有 
有 x 坐标 值 的 最 小 值 ， 则 a 是 种 群 最 高 等 级 中 的 个 体 。 提出 算法 在 非 支配 排序 时 ， 仅 处 理 当 前 种 群 中 最 高 等 级 个 
推理 2 若 存 在 个 体 5(x,y,z) 的 y 坐标 值 ， 是 当前 种 群 中 所 ” 体 ， 同 时 检测 是 否 选 足下 一 代 个 体 ， 当 发 现 选 够 时 ， 即 终止 算 
体 
前 


二 


有 yy 坐标 值 的 最 小 值 ， 则 5 是 种 群 最 高 等 级 中 的 个 体 。 法 ， 大 幅 减 少 了 运算 时 间 。 而 经 典 非 支配 解 排序 算法 需要 对 所 
推理 3 若 存 在 个 体 c(xy,z) 的 z 坐标 值 ， 是 当前 种 群 中 所 有 个 体 进 行 排序 ， 当 所 有 个 体 都 分 配 等 级 后 ， 再 选择 下 一 代 ， 
有 = 坐标 值 的 最 小 值 ， 则 是 种 群 最 高 等 级 中 的 个 体 。 所 以 用 时 较 多 。 第 一 代 NSGA 算法 的 时 间 复 杂 度 是 O(mm’), 虽 
推理 4 若 存 在 个 体 d(x,y,z) 到 原点 距离 的 值 ， 是 当前 种 群 ” 然 Deb 之 后 提出 了 改进 的 NSGAII, 但 排序 时 间 复 杂 度 是 OOxzz?)， 
中 任意 个 体 到 原点 距离 值 的 最 小 值 ， 则 4 是 种 群 最 高 等 级 中 的 。 原因 是 算法 在 运行 时 ， 要 对 所 有 个 体 分 配 等 级 后 ， 才 能 再 选 出 


个 体 。 下 一 代 个 体 ， 存 在 元 余 操 作 。 
证 明 著名 算法 ENSID 在 排序 时 ， 有 两 个 主要 步 又: 1) 把 种 群 按 


采用 反 证 法 ， 已 知 当前 种 群 ， 所 有 个 体 中 x 的 坐标 值 最 小 第 一 个 目标 函数 值 升序 排列 , 时 间 复 杂 度 为 O(n109”); 2) 对 每 个 
为 xnn， 如 果 个 体 a 坐标 是 aCxmin,aw,az), 则 当前 种 群 中 不 存在 一 解 分 配 等 级 。ENS-SS 在 非 支配 解 排序 时 ,最 坏 情况 是 种 群 中 所 


个 个 体 /， 使 得 7 的 横 坐 标 庆 小 于 xmin。 有 个 体 都 属于 同一 等 级 ， 此 时 时 间 复 杂 度 是 O(mn?”)， 最 好 情况 

本 文 假设 ， 存 在 一 个 个 体 71， 使 得 该 个 体 支配 个 体 是 种 群 个 体 属于 i 个 等 级 左右 ,时间 复杂 度 是 O(mn Vn )。 在 
aCxmimaywqaz)， 则 根据 Pareto 相关 原理 ， 可 得 到 [<xwmin, Ly<ay ENS-BS 进行 非 支配 解 排序 时 ,最 坏 情 况 和 ENS-SS 相同 , 时 间 
P<azs 明 显 ， 该 结论 与 已 知 条 件 矛 盾 ， 即 不 存在 个 体 1 支配 a， 复杂 度 是 O(mn”)， 最 好 的 情况 是 所 有 解 的 等 级 都 不 相同 ， 时 间 


则 得 出 个 体 a 属于 当前 最 高 等 级 个 体 , 证 毕 。 其 他 推理 2、3 和  ” 复杂 度 为 O(mnlog”)。 


4 的 证 明 过 程 与 此 类 似 。 由 此 可 知 ， 提 出 方法 中 首先 确定 4 个 本 文 提 出 的 算法 中 , 最 好 情况 是 第 一 等 级 个 体 数 量 > N。 等 
特殊 个 体 ， 一 定 是 当前 种 群 中 最 高 等 级 的 个 体 。 于 NN 时 ,直接 将 该 等 级 所 有 个 体 进入 下 一 代 ; 大 于 时 , 进行 拥挤 


距离 操作 ， 选 择 部 分 个 体 进入 下 一 代 ， 剩 余 所 有 个 体 不 再 进行 


和 1 八 、 | 
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等 级 分 配 ,所 需 时 间 复 杂 度 是 O(mn)， 算法 按 目 标 函 数值 排序 的 。” ”函数 描述 如 表 1 所 示 。FMODE 在 实验 中 参数 的 取 值 如 下 : 种 
时 间 复 杂 度 为 O(n1og”)， 计 算 到 原点 距离 时 间 复 杂 度 为 O(n)， 群 规模 N=200, 变 异 概率 =0.5, 交叉 概率 Cr=0.2, 由 于 不 同 函 
故 总 时 间 复 杂 度 为 O(mn)+2*O(nlog")+O(n)。 最 坏 情况 下 ,本文 ” 数 中 接近 理论 最 优 解 时 所 需 的 进化 代数 不 同 ， 设 置 不 同 的 进化 
算法 与 NSGAII 相同 ， 要 对 所 有 个 体 分 配 等 级 ， 时 间 复 杂 度 为 。” 代数 如 表 1 所 示 ， 选 变异 方式 为 “DE/Rand/1/bin”。 对 比 算法 为 
O(mm?), 总 时 间 复 杂 度 为 OGmn?)+2*O(nlog")+O(n), 在 进化 初始 著名 的 NSGAI 和 ENS 算法 , 他 们 的 参数 依据 算法 的 原文 献 进 
阶段 ， 初 始 种 群 个 体 分 布 较为 随机 ， 将 要 处 理 多 个 等 级 才能 i 行 设置 。 实验 环境 : 操作 系统 Win10, 处 理 器 Intel®) Core™i5- 
足 个 体 进入 下 一 代 。 但 是 ， 在 进化 中 后 期 ， 从 前 几 个 等 级 就 能 。” 4590CPU@ 3.3GHz ， 软 件 平台 Matlab R2015b 。 由 于 


选 足下 一 代 ， 从 而 提高 了 计算 效率 。 IGD(inverted generational distance)P3 包 含 了 多 样 性 和 收敛 性 , 且 
3 ”仿真 实验 能 体现 算法 性 能 ,所 以 本 文选 用 IGD 作为 测试 指标 , 详 见 式 (11) 
本 所 示 。APx 代 表 Pareto 理论 最 优 解 集 ，P 代 表 求 得 的 最 优 解 集 ，v 

3.1 实验 环境 及 参数 设置 e P，d(vwPD 代 表 v 和 P 之 间 的 最 小 欧 拉 距 离 。 

为 了 测试 算 ; 4 有 效 性 ， 采 用 国际 标准 测试 问题 

为 了 测试 算法 FMODE 的 有 效 性 ， 采 用 国际 标准 测试 问题 16D(P*.p)= Zsa) GD 
DTLZIL-DTLZ6， 这 些 问题 的 Pareto 前 沿 具 有 一 定 代 表 性 ， 有 具体 | | 

表 1 6 个 三 目标 测试 函数 
问题 ” 取 值 范围 标 函 数 维度 


DTLZ1 xe[0,1] {1 (=) (ts(%,)) m+4 


之 
g (x,) -oj + >》， (Gs -0.5) —cos(20x(% -03)) 


A(*)=(1+ s(%))eo 到 je 至 | sin (| 


DTLZ2 xel[0,1] m+9 


f(x) -rs(w))eo( 笃 jo[ 党 jeo [转世 筷 志 
f(D)=(1+ 8(% je 至 je 到 ) . oo 的 ja| 鸭 二 
) 


DTLZ3 xel[0,1] | m+9 


g (x,) -oj + >》， (Gs -0.5) —cos(20x(% 05) 


EXm 
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续 表 
问题 取 值 范围 标 函 数 维度 


DTLZ4 xe[0,1] m+9 


g(%,)= > (x -0.5) ， 2Q=100 


f(x) =(1 + 8(%, jos 至 je 至) co 5 到 
0 At 


f(x) IE : 


IE Ea 
/res)sin 经 ] 
加 -re 8) ) Vi (m+) 


fn 1(X)=%, 1 
DTLZ6 xe[0,1] fa(*)= J fs fi 8) m+9 


Ds 


Ni EXn 


8 (%,)= fs 


h(fisfssiis fg)=m a 


=l 


3.2 ”运行 速度 对 比 在 表 2 和 3 中 ，NP 为 种 群 个 数 , 下 为 特定 等 级 个 数 ， 括 弧 
为 验证 提出 排序 方法 的 高 效 性 , 本 节 给 出 与 NSGAIT、 ENS ” 外 是 排序 耗 时 ， 单 位 为 s， 括 弧 内 是 比较 次 数 ， 表 2 使 用 具 
非 支 配 解 排序 方法 在 运行 速度 和 比较 次 数 上 对 比 , 结果 如 表 2、 ”特定 等 级 个 数 的 种 群 , 并 设置 等 级 个 数 为 YNP 左右 , 在 此 等 级 
3 所 示 。 表 中 时 间 单 位 是 s， 比 较 次 数 单位 是 次 ， 在 不 同 种 群 规 。 下 ENS-SS 表现 出 最 优 性 能 , 详 见 文献 [11]。 
模 下 分 别 测试 ， 测 试问 题 都 是 三 目标 优化 问题 。 算 法 ENS_SS 通过 表 2 和 3 可 以 看 出 ， 本 文 的 排序 算法 在 个 体 数量 较 少 
和 ENS_BS 分 别 代 表 的 是 采用 DE 进化 框架 结合 ENS_SS 排序 ”时 排序 耗 时 与 其 他 三 种 算法 相差 不 大 。 当 种 群 个 体 较 多 时 ， 本 
和 ENS_BS 排序 。 需 要 说 明 的 是 ， 在 本 节 测 试 中 ， 四 种 算法 对 文 算法 排序 耗 时 略 多 于 ENS 方法 , 但 远 小 于 NSGAII 的 排序 耗 
种 群 中 所 有 个 体 仅 计算 等 级 操作 ， 不 做 拥挤 距离 计算 。 即 使 对 ”时 。 简 言 之 ,FMODE 表现 计算 效率 略 进 于 ENS 算法 ， 但 都 显 
所 有 个 体 都 计算 等 级 处 理 ， 本 文 算法 每 次 得 到 当前 最 高 等 级 个 。” 著 优 于 NSGAII 算法 。 由 于 ENS 是 近年 来 提出 的 著名 算法 , 与 
体 ， 在 后 续 等 级 分 配 中 ， 则 不 需要 再 和 已 分 配 过 等 级 的 个 体 对 ”其 对 比 实验 结果 显示 出 其 的 优越 性 ， 也 说 明了 提出 算法 在 计算 
比 ， 理 论 上 与 ENS 相差 不 大 ， 但 比 NSGAII 效率 更 高 。 速度 上 的 高 效 性 。 
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表 2 四 种 排序 方法 对 于 具有 NP 等 级 种 群 的 对 比 测试 :排序 耗 时 秒 (比较 次 数 ) 
NP,F 
法 
100, 10 500,22 1000,32 2000,45 5000,71 10000, 100 
0.064292 1.2045 4.7786 19.7604 143.8898 713.4332 
NSGAII 
(4950) (124750) (499500) (1999000) (12497500) (49995000) 
0.0011174 0.0055482 0.020251 0.041009 0.15368 0.41325 
ENS SS 
(900) (10972) (30780) (88080) (350050) (990000) 
0.00090068 0.003208 0.011526 0.02105 0.079892 0.18858 
ENS BS 
(646) (6999) (18415) (51171) (195664) (543093) 
0.020339 0.028142 0.079487 0.14851 0.61336 1.7341 
FMODE 
(879) (10921) (30704) (87975) (349889) (989758) 
表 3 四 种 排序 方法 对 于 随机 生成 种 群 的 对 比 测试 :排序 耗 时 秒 (比较 次 数 ) 
NP 
法 
100 500 1000 2000 5000 10000 
0.012643 0.30389 1.2565 5.1451 33.629 152.3857 
NSGAII 
(4950) (124750) (499500) (1999000) (12497500) (49995000) 
0.0007792 0.0068237 0.015762 0.044478 0.20465 0.56412 
ENS SS 
(1237) (16263) (49164) (147723) (627683) (1910611) 
0.00079018 0.0078485 0.021014 0.064339 0.32262 1.1166 
ENS BS 
(1451) (22916) (76597) (250421) (1179679) (4015154) 
0.012497 0.034995 0.072663 0.22605 0.96877 2.9161 
FMODE 
(1395) (16994) (51197) (151049) (636641) (1927848) 
把 ( 洒 性 能 天 Sr a 
3.3 总 体 性 能 对 比 表 4 四 种 算法 求解 最 小 化 三 目标 问题 
本 节 给 出 FMODE 与 经 典 算 法 总 体 性 能 比较 ， 对 比 结果 如 DTLZ1-DTLZ6 上 IGD 均值 测试 对 比 
表 4、5 所 示 ， 算 法 性 能 评价 采用 IGD 指标 ， 指 标 详细 说 明 见 ee 
文献 [16]， 表 中 的 最 优 结果 ] 粗 体 显示 , 测试 结果 为 五 次 独立 运 函数 ENS _ SS00 ENS BS NSGAIID] 本 章 算法 (EMODE) 
行 的 才 司 上 时 给 由 月 小 在 六 个 测 汪 卫 当 和 运 管 时 下 
行 的 均值 。 同 时 给 出 所 有 算法 在 六 个 测试 函数 上 的 运算 时 间 和 TO TD on TD TO 
线 ， 结 果 如 图 2~7 所 
进化 曲线 ， 结 果 如 图 2~7 所 示 。 DTLZ2 4.44E-02 4.45E-02 4.37E-02 4.30E-02 
习 区 2~ 向 疆 里 > | 汗 丰 求 六 个 
从 表 4、5 和 图 2~7 的 结果 可 以 看 出 ， 四 种 算法 在 求 六 1 DTLZ3 ”4.49E-02 4.46E-02 4.44E-02 4.45E-02 
函数 时 进化 曲线 大 致 相似 ， 其 中 ENS-BS 和 ENS-SS 两 种 算法 Driz4 439E-02 445E-02 4.52E-02 4.54E-02 
ER A a 
稍微 占 优 ， 这 得 益 于 ENS 算法 把 见 余 计算 尽量 降低 。 在 求 角 DTLZ5 2.70E-03 2.57E-03 2.97E-03 2.96E-03 
ya : 用 时 稍 旦 
DTIZ1 和 DTLZ3 函数 时 , FMODE 进化 用 时 稍 多 于 ENS, 但 是 DTLZ6 432E-02 4.02E-02 8.43E-02 7.71E-02 
远 少 于 NSGAII。 在 求解 DTLZ4 函数 表现 上 ，FMODE 性 能 优 
上 寸 稍 筑 人 六 Y WN Ne 
于 NSGAH， 所 用 耗 时 稍微 差 于 NSAGII。 在 解决 DTLZ2 和 表 5 四 种 算法 求解 最 小 化 三 目标 问题 
DTLZ6 时 ，FMODE 表现 性 能 与 对 比 算法 基本 一 致 ， 但 用 时 稍 DTLZ1-DTLZ6 进化 完成 所 用 时 间 均 值 测试 对 比 
差 与 NSGAII, 特 别 在 解决 DTLZ6 时 ,FMODE 用 时 明显 较 多 ， i ea ne 
JR 入 AN > 左 求 解 沙 - pp 个 - 取 侍 在 季 困 炎 
深入 分 析 原 因 ， 发 现在 求解 该 问题 时 ， 大 部 分 | 体 都 聚集 在 第 ENS_SS0D ENS BS NSGAIIDI 本 章 算 法 (FMODE) 
一 人 了 芒 引 S 第 一 笠 《 A » 慎 4 二 全 
等 级 ， 所 以 FMODE 需要 对 第 “等 级 的 所 有 个 体 进行 计算 处 DTLZ1 2.09E+01 2.13E+01 3.36E+01 2.56E+01 
ee es 
理 ， 无 法 体现 其 在 计算 的 同时 选择 下 一 代 个 体 的 优势 。 算 法 计 DTLZ2 2.91E+01 2.71E+01 2.69E+01 2.65E+01 
得 到 的 前 端 与 理论 前 端 对 比 见 图 
算得 到 的 前 端 与 理论 前 端 对 比 见 图 8 所 示 。 DTLZ3 4.94E+01 5.55E+01 8.64E+01 6.11E+01 
此 ES 亚 价 起 村 柱子 日 
总 体 而 言 , 在 IGD 评价 指标 性 能 方面 ,提出 的 FMODE 算 DTLZ4 1.33E+01 1.30E+01 1.79E+01 1.75E+01 
2 间 二 上 日立 、 性 人 川 其 百 | 
法 的 与 ENS、NSGAII 相差 不 大 ， 人 性 能 表现 基本 相同 。 但 是 在 DTLZ5 2.01E+01 2.01E+01 1.98E+01 1.80E+01 
和 S 六 证 | 当 阳 1 关于 [= 仿古 
时 间 对 比方 面 ， 稍 微 差 于 ENS， 明 显 优 于 NSGAII， 从 而 验证 DTLZ6 2.15E+02 2.06E+02 2.30E+02 2.51E+02 
了 提出 算法 的 有 效 性 和 正确 性 . 
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表 6 四 种 算法 求解 最 小 化 三 目标 问题 
DTLZ1-DTLZ6 进化 完成 所 需 代数 测试 对 比 


DTLZ1 DTLZ2 DTLZ3 DTLZ4 DTLZS DTLZ6 


进化 代数 100 80 300 55 70 1000 


102 


101 


Mean Error (IGD) 
己 


二 
马 


102 
20 40 60 80 100 
FES 


Mean Error (TIME) 
5 


20 40 60 80 100 
FES 


(b) 
图 2 DTLZ1 进化 曲线 (a) 与 进化 耗 时 (b) 比 较 


Mean Error (IGD) 


Mean Error (TIME) 


10 20 30 40 50 60 70 80 
FES 


(b) 


图 3 DTLZ2 进化 曲线 (a) 与 进化 耗 时 (b) 比 较 


Mean Error (IGD) 


Mean Error (TIME) 
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(a) 
10< 
101 
10° 
(b) 
图 4 DTLZ3 进化 曲线 (a) 与 进化 耗 时 (b) 比 较 
10° 
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(b) 


图 5 DTLZ4 进化 
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化 耗 时 (b) 比 较 


201805.00399vVv1 


chinaXiv 


录用 稿 


Mean Error (IGD) 
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(a) 
图 6 DTLZ5 进化 
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E] 


Mean Error (TIME) 
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(b) 


线 (a) 与 进化 耗 时 (b) 比 较 
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线 (a) 与 进化 耗 时 (b) 比 较 
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图 7 DTLZ6 进化 
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图 8 DTLZ1-DTLZ6 理论 前 端 与 FMODE 算法 得 到 的 前 端 对 比 ， 红 
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色 为 理论 前 端 ， 黑 色 为 算法 计算 得 到 前 端 ，(a)~( 了 分 别 为 DTLZ1-DTLZ6 


3.4 算法 FMODE 的 参数 设置 
本 节 研 究 提出 算法 FMODE 的 参数 设置 ， 给 出 在 不 同 交叉 
概率 Cr、 变 异 因 子 F、 变 异 策略 的 组 合 对 算法 性 能 的 影响 。 


在 变异 方式 “DE/Rand/1/bin” 中 ， 变异 因子 F=0.5， 群体 规 
模 NP=200， 最 大 进化 代数 FES=1000 代 时 ， 对 交叉 概率 Cr 进 
行 不 同 的 取 值 测试 ， 其 IGD 均值 结果 如 表 7 所 示 。 通 过 表 7 可 
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以 发 现 ,“DE/Rand/l/bin” 变 异 方式 下 ， 缩 放 因子 F=0.5， 交 叉 
概率 Cr=0.2 或 0.3 时 效果 最 好 , 由 此 得 到 F=-0.5 和 Cr=0.2 作 为 ” [9] 
经 验 值 。 在 其 他 变异 策略 下 ， 不 同 的 交叉 和 缩放 因子 组 合 时 ， 
所 表现 的 性 能 都 差 于 “DE/Rand/1/bin” 变 异 ， 为 节省 版 面 其 数 
据 不 再 给 出 。 最 终 本 文 得 出 ， 提 出 算法 推荐 的 参数 设置 为 : [7] 
“DE/Rand/1/bin ”变异 ， 缩 放 因 子 F=0 .5， 交 叉 概 率 Cr=0.2。 


26 (8): 997 - 1003. 


Cui Xunxue, Li Qin, Tao Qing. Genetic algorithm for Pareto optimum-based 


route selection [J]. Journal of Systems Engineering and Electronics, 2007, 2 


(18): 360-368. 


王 勇 , 蔡 自 兴 ， 


周 育 人 , 等 . 约束 进化 算法 [加 . 软件 学 报 , 2009, 20 (1): 


11-29. 


表 7 DE/rand/1/bin”* 变 异 时 不 同 参数 值 的 IGD 均值 对 比 [8] 郑 建国 , 王 翔 , 刘 荣 辉 . 求解 约束 优化 问题 的 e -DE 算法 [J]. 软件 学 
缩放 因子 和 ae 报 , 2012, 23 (9): 2374-2387. 
交叉 概率 DTLZ1 DM i [9] Tang Suqin, Cai Zixing, Zheng Jinhua. A fast method of constructing the 
SEE 人 non-dominated set: arena’s principle [C]/ Proc of the 4th International 
=0.5,Cr=0. . -| . - 和 - 了 -| % -| 交 - 
Conference on Natural Computation. 2008: 391--395. 
F=0.5,Cr=0.2 .46E-02 3.92E-02 4.01E-02 4.01E-02 2.22E-03 3.45E-02 
BO Dd “Ado ods Woo oid dn .Cel [10] Mc Clymont K, Keedwell E. Deductive sort and climbing sort: new methods 
= 和 CT. 到 - -| - - « -| -| 
ee SO. den “dd deo. odd a for non-dominated sorting [J]. Evolutionary Computation, 2012, 20 (1): 1— 
P=0.5.Cr=0.5 .601E-02 4.52E-02 1.1SE+01 4.69E-02 2.06E-03 2.60E-01 26. 
F=0.5,Cr=0.6 1.72E-02 4.81E-02 4.86E+01 4.97E-02 2.17E-03 4.83E-01 U1] Ahang Xingyl, Tian Ye Cheng Ran, etal. An.efficlent approach to non: 
Re sy BO dD Ooo Voidy. Ws. dnl dominated sorting for evolutionary multi-objective optimization [J]. IEEE 
F=0.5,Cr=0.8 26E+00 5.19E-02 1.18E+02 5.55E-02 2.48E-03 1.55E+00 nn 
=0.5,Cr=0. -| ， “ 由 
[12] Xu Yulong，Fang Jian’an, Zhu Wu, et al. Differential evolution using a 
superior-inferior crossover Scheme [J]. Computational Optimization and 
4 ”结束 语 加 
Applications, 2015, 61 (4): 243-274. 
本 文 介绍 一 种 高 效 标 非 支配 解 排序 方法 ， 该 方法 能 降 [13] Das S$, Mullick S S, Suganthan P N. Recent advances in differential 


低 排序 操作 处 理 个 体 的 数量 , 且 在 排序 操作 的 同时 ,可 进行 选择 
下 一 代 个 体 的 操作 ， 较 大 提升 原始 排序 分 配 等 级 的 计算 效率 。 
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